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We introduce a generalization of Wegner's n-orbital model 
for the description of randomly disordered systems by replac- 
ing his ensemble of Gaussian random matrices by an ensemble 
of randomly rotated matrices. We calculate the one- and two- 
particle Green's functions and the conductivity exactly in the 
limit n — > oo. Our solution solves the CPA-equation of the 
(n = l)-Anderson model for arbitrarily distributed disorder. 
We show how the Lloyd model is included in our model. 

PACS numbers: 71-10.+X, 05.60.+W, 71.55.-i 



The treatment of physical systems with disorder 
presents one of the great challenges in statistical physics. 
Realistic multi-site models, like the Anderson model |Q, 
are in general unsolvablc. Exact calculations are only 
possible in one pathological special case - namely for the 
Lloyd model M with Cauchy-distributed disorder. Oth- 
erwise one has to use approximation techniques which re- 
duce the multi-site to single-site models and treat them 
self-consistently. The most prominent of these methods 
is the so-called coherent potential approximation (CPA), 
see, e.g., [|L 

Wegner] 4| introduced a generalization of the Anderson 
model by putting n electronic states at each site and 
describing the disorder by Gaussian random matrices in 
these electronic states. Whereas for n = 1 this reduces 
to the usual (unsolvable) Anderson model, the opposite 
limit n — > oo becomes exactly solvable. Interestingly, the 
solution of this multi-site model coincides with a special 
CPA-solution. 

In this letter we introduce a generalization of Wegner's 
model by replacing his ensemble of Gaussian random ma- 
trices by a more general ensemble of random matrices - 
and thus allowing arbitrarily distributed disorder. Nev- 
ertheless - by using recent results of Voiculescu and 
Speicher |7]|| on the mathematical concept of "freeness" 
- we are still able to calculate the Green's functions of this 
model in the limit n — + oo exactly. Again our solution 
solves the CPA-equation of the (n = l)-Anderson model. 
Hence our model can be regarded as a rigorous mean-field 
model for CPA for arbitrarily distributed disorder. 

Let us first recall Wegner's model. He considers a d- 
dimensional lattice where at each site r there are n elec- 
tronic levels \ra) numbered by a = 1, . . . , n. The inter- 
action is governed by a Hamiltonian of the form 

H = H + H U (1) 

where Hq is deterministic, translational- invariant, and 
diagonal in the electronic levels, 



H = ^ v r -r'\ra)(r'a\, 

and Hi describes the site-diagonal disorder 

Hi= £ -±=ff\ra){r!3\. 
a v n 



(2) 



(3) 



There, f r = ((l/v 7 ^)/"^)™ /3 _ 1 are Gaussian random 
matrices and the entries of f r and of f r i are indepen- 
dent for r ^ r' . (This is the site-diagonal model of Weg- 
ner, we will not treat his "local gauge invariant" model 
here.) By using techniques for calculating moments of 
Gaussian random matrices ||[l0]], Wegner was able to 
calculate the one- and two-particle Green's function in 
the limit n — > oo. 

To explain our generalization let us diagonalize the 
random matrices f r in the form f r — u r f r u*. Thus the 
ensemble of the f r is determined by an ensemble of diag- 
onal matrices f r and an ensemble of unitary matrices u r . 
In the limit n — > oo, Wegner's original formulation is re- 
covered if the diagonal f r are taken as deterministic ma- 
trices having Wigner's semi-circle || as eigenvalue distri- 
bution and if the u r are random unitary matrices, given 
by the canonical invariant or Haar measure on U(n), such 
that u r and u r > are chosen independently from U(n) for 
r ^ r'. Since this ensemble is invariant under indepen- 
dent rotations at different sites we may also replace the 
different f r by one single (not necessarily diagonal) /, i.e. 
we have f r = u r fu* with the ensemble of tiy as stated 
above. 

The advantage of this reformulation of Wegner's model 
is that now a generalization is obvious: We are a priori 
totally free in the choice for /. Thus our model is given 
by the following Hamiltonian 



H = 



E 



\ra)(r'c 



E ( u rfK)a,!3\ra)(rf3\, 

r,a,0 



(4) 



where / is a deterministic hermitian n x n-matrix and 
u r G U(n) are random unitary matrices, chosen indepen- 
dently for different sites. This means that we act at each 
site r with a copy f r := u r fu* of the given operator /, 
but that the basis for f r and the basis for f r i are rotated 
randomly against each other for all pairs of different sites 

The possibility for treating this model in the limit n — > 
oo arises from the important observation of Voiculescu j^] 
(see also pT[) that there is a connection with his concept 
of free random variables Jl3|,|6|Jl^] : Denote by (. . .) cns the 



1 



average over our ensemble of random unitary matrices 
and let 



1 



(al- 



count the averaged eigenvalue distribution of our n x 
n-matrices. Then we have for all m S N and for all 
polynomials pi,... ,p m with (pi(f)) = (i = 1, . . . , m) 
in the limit n — > oo that also 



<Pl(/r(l))p 2 (/r( 2 )) ■ ■ ■Pm(/r(m))) = 



(6) 



for all sequences of indices r(l), . . . , r(m), where all con- 
secutive indices are different, i.e. where r(i) ^ r(i + 1) 
for all i = 1, . . . , m — 1. An example of ([|) is 



(Pi(/i)P2(/2)P3(/i)P4(/ 2 )) = 0, 



(7) 



whereas for sequences with coinciding neighbouring in- 
dices one gets non-vanishing results like 

(pi(/l)p 2 (/2)p 3 (/2)M/l)} = (8) 
(pi(/l)P4(/l))<P2(/ 2 )p3(/ 2 )) = (pl(/)P4(/)><P2(/)p 3 (/)>. 

The property (^) - in the mathematical literature recently 
introduced by Voiculescu fll3|,|6| under the name of "free- 
ness" - allows, as in Wegner's case, to calculate all mixed 
moments of the matrices f r (like, e.g., in (||)) and thus to 
derive exact expressions for the Green's functions of our 
model in the limit n — > oo. Whereas usually one gets an 
infinite hierarchy of equations for averaged quantities, it 
is exactly the property of freeness (^|) which closes our 
equations. For an effective handling of the calculations 
one needs the i?-transform machinery of Voiculescu fbipl 
and the concept of non-crossing cumulants of Speicher 
j^D- The concrete calculations will be published else- 
where, here we only want to give the results. 

The most important quantity is (the diagonal part of) 
the one-particle Green's function (IPG) 



G(z) 



n 
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H 



ra 



(9) 



whose spectral function yields the density of states. Note 
that, due to the translation- invariance of our Hamilto- 
nian, G{z) is independent of r and hence the local and 
the global density of states coincide. If we denote sim- 
ilarily by Go and G\ the IPG of the deterministic and 
the disorder part of H, respectively, i.e. 



(10) 




(11) 



then the crucial property (Q) allows us to derive the fol- 
lowing self-consistent equation for G(z): 



G(z) = G (z - i?i(G(z))), 



(12) 



where R\(z) is a kind of self-consistent self-energy for H\, 
namely it is defined by 



Gx{z) = 



1 



z-R x {G l {z)) 



(13) 



with i?i(0) = 0. In the special case of a semi-circle eigen- 
value distribution for /, G\ is given by 



Gi{z) = 



z-Vz 2 



4ct 2 
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z-a 2 G x {z) 



(14) 



showing that R\{z) = a 2 z. In this case our formula (|l 
reduces to Wegner's result P,|l5| 



G(z) = G (z-ct 2 G(z)). 



(15) 



Given G(z) and Ri(z), one can derive linear self- 
consistent equations for the one-particle and two-particle 
Green's functions. Again, the freeness property @ is 
the essential ingredient for closing the equations. For 
the Fourier-transform of the one-particle Green's func- 
tion one obtains 



G(q;z) 



1 



z-v(q)-R 1 (G(z)) 



(16) 



with v r = J q v{q)e lqr where J q = J 1BZ d d q, V being 
the volume of the first Brillouin zone (1BZ). For R\(z) = 
a 2 z this reduces again to Wegner's result |B,^5[- 

The two-particle Green's function determines the con- 
ductivity as a velocity-velocity correlation function. Al- 
though in general the two-particle Green's function con- 
tains an additional term, one obtains by the symmetry 
argument u_ r = v r that only the product of the one- 
particle Green's functions contribute. Thus one finds for 
zero temperature and in the dc-limit uj — > 



a(w,0) 



2ne 2 B 
V 



where 
B := 



g 2 {R 1 {G);E F ) 



\Vv{q)[< 



(17) 



q \v(q) - ((E F +lu + iQ+)\ 2 \v(q) - C(Ef + *0+)| 2 
with C(z) = z — R\(G{z)) and the spectral function 



q{R 1 {G);E) = — lmR 1 (G(z = E + iO + )). 



(18) 



For R\(z) = a 2 z this reduces to Wegner's result in a form 
given by Khorunzhy and Pastur [|15| . 
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The generality of our solution may be seen from the 
fact that it includes also the Lloyd model as a special 
case. Namely, choose a Cauchy-distribution with param- 
eter 7 as eigenvalue distribution for /. Then 



1 r + °° 

Gx{z) = - 

T J-oo 



7 



1 



Y + t A z - t 



-dt = 



1 



z — ij 



(19) 



hence R\(z) = vy. This shows that the Cauchy- 
distribution behaves in all relevant aspects like an imagi- 
nary (^-distribution. For ^-distributions, however, there is 
no difference between our model and the original Ander- 
son Hamiltonian (see also [[l6 17 1). Since the Anderson 
model with Cauchy-distributed disorder is nothing else 
but the Lloyd model, the latter is included in our inves- 
tigations. In particular, we recover from ( |l~2"l) the IPG of 
the Lloyd model § 



G(z) = G (z-ij). 



(20) 



One surprising feature of the Wegner model is that its 
solution coincides with a special CPA-solution. This gen- 
eralizes also to our model. In general, the CPA-solution 
for the Anderson model with single-site random variable 
X is given by the two equations 



G{z)=G {z-H{z)) 



and 



X - £(z) 
l-(A-£(z))G(z) 



= 0. 



(21) 



(22) 



The first of these equations coincides with our solution 
(H) if we identify £(z) = R x (G(z)). By using the equiv- 
alent form of (|l]), namely 



Gx{Ri{z) + z- 1 ) 



(23) 



it can be checked that our solution fulfills also (|22j) if we 
choose as distribution for X the eigenvalue distribution 
of /. Note that we specify our model rigorously in the 
beginning and that we are able to calculate all quantities 
without any further approximation. Thus our multi-site 
model is a rigorous mean-field model for the usual single- 
site CPA. One should also note that our previous remarks 
about the Lloyd model can now be taken as an explana- 
tion for the well-known fact that CPA is exact for the 
Lloyd model. 

As an instructive example of our formalism let us con- 
sider the one-dimensional lattice with nearest-neighbour 
interaction and binary site-diagonal disorder for the spe- 
cial case v = a = 1, where 2v is the half-band- width 
and a 2 is the variance of the disorder. Then we have 
G (z) = l/Vz 2 - 1 and R x {z) = (Vl + z 2 -l)/2z, which 
yields as a solution of (JT3) 



G(z) = 



8 

6Vz 2 -2 + 2z' 



(24) 



This provides 

g(Ri (G);E) = -^-\/2- E 2 Q(V2-E), (25) 

i.e. together with (|l~7| ) a finite conductivity everywhere 
inside the band [—V%, V2]- 

Since we can prescribe in our model an arbitrary dis- 
tribution for / - or for X in the CPA- formulation -, it 
is clear that our formalism is capable of covering a lot 
of quite different examples, like, e.g., a superposition of 
binary noise or, more general, the recently investigated q- 
noise p"8[ |. Furthermore, our description using the theory 
of free random variables and the notion of non-crossing 
cumulants allows a straightforward generalization to the 
case of dynamical disorder and thus promises to give a 
rigorous model for dynamical CPA. These subjects will 
be pursued further in forthcoming investigations. 

We acknowledge helpful discussions with Franz Weg- 
ner and Petr Chvosta. This work was supported by the 
Deutsche Forschungsgemeinschaft (R.S.). 
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